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ABSTRACT 

Let  r  be  a  collection  of  n  (possibly  intersecting)  "red"  Jordan  arcs  of  some 
simple  shape  in  the  plane  and  let  T'  be  a  similar  collection  of  m  "blue"  arcs. 
We  present  several  efficient  algorithms  for  detecting  an  intersection  between 
an  arc  of  T  and  an  arc  of  T'.  (i)  If  the  arcs  of  V  form  the  boundary  of  a 
simply  connected  region,  then  we  can  detect  a  "red-blue"  intersection  in  time 
0  (Xs{''TT')  log^  m  +  {Xs{m)  +  n)  log(n  +  m)  j  where  \s{m)  is  the  (almost-linear) 
maximum  length  of  (m,  5)  Davenport-Schinzel  sequences,  and  where  s  is  a  fixed 
parameter,  depending  on  the  shape  of  the  given  arcs.  Another  case  where  we  can 
detect  an  intersection  in  close  to  linear  time  is  when  the  union  of  the  arcs  of  T 
and  the  union  of  the  arcs  of  F'  are  both  connected,  (ii)  In  the  most  general  case, 
we  can  detect  an  intersection  in  time  O  f  (m-«/Aj(n)  -I-  7i>/Aj(m))  log^'^(m  +  n)) . 
For  several  special  but  useful  cases,  in  which  many  faces  in  the  arrangements  of 
F  and  F'  can  be  computed  efficiently,  we  obtain  randomized  algorithms  which 
are  better  than  the  general  algorithm.  In  particular  when  cdl  arcs  in  F  and 
F'  are  line  segments,  we  obtain  a  randomized  0{{m  +  n)^^^'^^)  intersection  de- 
tection algorithm,  for  any  e  >  0.  We  apply  the  algorithm  in  (i)  to  obtain  an 
0(A,(n)log^  n)  algorithm  (for  some  smcdl  s  >  0)  for  planning  the  motion  of  an 
n-sided  simple  polygon  around  a  right-angle  corner  in  a  corridor,  improving  a 
previous  0{n^)  algorithm  of  [MY86],  and  to  derive  an  efficient  technique  for  fast 
collision  detection  for  a  simple  polygon  moving  (translating  and  rotating)  in  the 
plane  along  a  prescribed  path. 


Work  on  this  paper  has  been  supported  by  Office  of  NavaJ  Research  Grant  N00014-82- 
K-0381,  by  National  Science  Foundation  Grant  No.  NSF-DCR-83-20085,  by  grants  from  the 
Digital  Equipment  Corporation,  and  the  IBM  Corporation.  Work  by  the  second  author  has 
also  been  supported  by  a  research  grant  from  the  NCRD  -  the  Israeli  National  Council  for 
Research  and  Development. 


1      Introduction 

Consider  the  following  problem.  Let  F  be  a  collection  of  n  (possibly  intersecting) 
"red"  Jordan  axes  in  the  plane,  and  let  T'  he  a  similar  collection  of  m  "blue" 
arcs.  Does  any  red  arc  intersect  any  blue  arc?  Suppose  the  arcs  have  relatively 
simple  shape;  e.g.  suppose  they  axe  line  segments,  or  x-monotone  algebraic  arcs  of 
some  small  fixed  degree,  so  that  in  particular  each  pair  of  them  intersects  in  at 
most  some  fixed  number  5  of  points.  There  is  a  simple  way  to  detect  a  red-blue 
intersection,  using  a  standard  sweep-line  algorithm,  such  as  that  of  Bentley  and 
Ottmcuin  [B079].  However,  this  algorithm  runs  in  time  0{{m  +n  +  t)log(m  +  n)), 
where  t  is  the  number  of  red-red  axid  blue-blue  intersections  that  the  algorithm  has 
to  sweep  through  before  detecting  a  red-blue  intersection.  Since  t  can  be  quadratic 
(i.e.  ^{mn))  in  the  worst  case,  the  algorithm  is  not  generally  efficient,  and  in  fact 
might  probably  be  inferior  to  the  naive  0{mn)  algorithm  that  checks  all  possible 
pairs  of  a  red  axe  axid  a  blue  axe  for  intersection. 

In  this  paper,  we  present  several  efficient  algorithms  for  detecting  a  red-blue 
intersection,  which  avoid  the  overhead  of  having  to  examine  many  red-red  or  blue- 
blue  intersections.  The  easiest  case  is  when  the  red  axes  in  T  do  not  cross  one  axiother 
and  form  the  boundary  of  a  simply  connected  region  (e.g.  segments  bounding  a 
simple  polygon).  In  this  ease,  we  can  detect  a  red-blue  intersection  in  time 

O  (A,+2(m)  log^  m  -I-  {\,+2{rn)  +  n)  log(m  +  n))  , 

where  Aj(m)  is  the  maximum  length  of  an  (m,  s)  Davenport-Schinzel  sequence  (see 
[HS86],  [ASS87]),  which  is  almost  linear  in  m  for  any  fixed  value  of  s;  specifically 

Ai(m)  =  m,      A2(m)  =  2m  —  1     (trivial) 

A3(m)  =  e{ma{m))     [HS86] 

\^(m)  =  0(m-2'*<'"))     [ASS87] 

A2,+2(m)  =  0{m  ■  2-'('")(i+°(i)))  for^  >  1     [ASS87] 

A2.+3(m)  =  0(m  •  2"'('")'°8°<"'^(^+°(i)))  for5  >  1     [ASS87] 

A2,+2(m)  =  fi(m-2"('*'('"»)for5>  1     [ASS87] 

where  a{m)  is  the  extremely  slowly  growing  functional  inverse  of  AckermEuin's  func- 
tion. Thus  we  obtain  a  close-to-lineax  solution  to  the  intersection  detection  problem 
in  this  special  case.  (In  the  more  special  case,  where  each  red  arc  and  each  blue  circ 
is  a  straight  line  segment,  we  can  use  the  ray-shooting  technique  of  [CG85]  to  detect 
an  intersection  of  each  red  segment  with  the  blue  polygon  in  O(logn)  time,  after 
0(n  log  log  n)  preprocessing,  thus  obtaining  an  improved  eJgoritlun  for  intersection 


detection.   However,  for  general  arcs,  no  such  efficient  "axc-shooting"  procedure  is 
known.) 

Our  aJgorithm  uses  a  recent  result  of  [GSS88]  which  shows  that  the  combinatorial 
complexity  of  a  single  face  in  an  arrzingement  of  n  axes,  each  pair  intersecting  in 
at  most  s  points,  is  0{\3+2{n)),  eind  that  such  a  face  caii  be  computed  in  time 
0(A,+2(")log^ '^)-  We  also  extend  our  algorithm  to  arbitrary  collections  F,  T'  of 
intersecting  axes,  provided  the  union  of  all  red  axes  is  a  connected  set,  as  is  the 
union  of  all  blue  axes;  in  this  extension  the  algorithm  runs  in  time  0{\t+2{m  + 
n)log^(m  +  n)). 

These  results  have  several  applications  in  motion  plaxming  and  collision  detection 
for  a  simple  polygon.  For  exaxnple,  Maddila  and  Yap  [MY86]  study  the  problem 
of  moving  (traxislating  and  rotating)  2m  n-sided  simple  polygon  B  airound  a  right- 
angle  comer  in  a  corridor.  They  show  that  the  problem  can  be  reduced  to  testing  a 
single  canonical  motion  of  B  for  intersection  with  the  corridor  walls,  which  in  turn 
can  be  reduced,  by  studying  the  problem  in  a  coordinate  fraxne  attached  to  B,  to 
testing  whether  B  intersects  some  resulting  collection  of  algebraic  axes  (describing 
the  relative  motion  of  the  corners  and  walls  of  the  corridor).  Our  intersection 
detection  algorithm  cam  then  be  used  to  obtain  an  0(A,(n)log^  n)  motion  planning 
algorithm  (for  some  smeill  s  >  0),  considerably  improving  the  Oin"^)  algorithm  given 
in  [MY86].  More  generally,  we  can  use  our  procedure  to  obtain  fast  algorithms  for 
testing  a  prescribed  motion  of  a  simple  polygon  (or  "curved"  polygon)  B  for  collision 
with  a  given  collection  of  obstacles.  Under  reasonable  assumptions  on  the  simplicity 
of  the  motion  of  B  (e.g.  along  algebraic  paths  of  low  degree)  and  of  the  shape  of 
the  boundaxies  of  the  obstacles  and  of  B,  we  can  obtain  close  to  lineax  collision 
detection  procedures. 

In  the  general  case  of  the  red-blue  intersection  detection  problem,  our  algorithms 
axe  less  efficient,  but  axe  still  significantly  better  than  quadratic.  In  a  nutshell, 
our  first  algorithm  is  efficient  because  it  can  restrict  the  problem  of  detecting  a 
red-blue  intersection  over  a  single  face  in  the  blue  or  the  red  axrangement.  For 
arbitraxy  red  and  blue  arrangements  this  is  not  possible  in  genered,  and  one  needs 
to  calculate  and  seaxch  for  an  intersection  in  many  such  faces.  We  then  face  the 
problem  of  obtaining  sheirp  bounds  for  the  complexity  of  many  such  faces,  and  of 
their  efficient  calculation.  For  general  arcs,  when  such  procedures  axe  not  avail- 
able, we  present  a  deterministic  algorithm  which  detects  a  red-blue  intersection 
in  time  O  {^{m^Xs+'iin)  +  n^J X,+2im))  log^-^(m  -|-  n)).  For  certain  special  types  of 
axes,  for  which  many  faces  in  their  axrangements  can  be  computed  efficiently,  we 
obtain  improved  (albeit  randomized)  algorithms.  For  example,  if  all  arcs  axe  line 
segments  or  unit  circles,  then  we  can  detect  an  intersection,  with  high  probabihty, 


in  time  0{{m  +  n)''/3+')  for  any  e  >  0,  exploiting  the  results  of  [CEGSW88]  [EGS88] 
[EGH*88]. 

We  believe  our  algorithms  can  be  extended  to  obtain  all  red-blue  intersections  in 
an  output-sensitive  fashion;  if  all  arcs  are  line  segments,  then  we  have  an  algorithm 
which  produces  all  k  red-blue  intersections  in  time  O  {{rn^/n  -\-  ny/m  +  k)  log(m  -|-  n)), 
using  the  ray-shooting  technique  of  [CG85]. 

The  paper  is  organized  as  follows.  Section  2  presents  the  efficient  algorithm  for 
the  case  of  a  simple  (curved)  polygon,  and  Section  3  gives  applications  of  this  algo- 
rithm to  motion  planning  and  collision  detection.  Section  4  describes  algorithms  for 
general  red-blue  intersection  detection,  and  mentions  a  few  additional  applications. 
Finally  in  Section  5  we  mention  some  open  problems. 

2      Intersection  between  a  simple  polygon  and  arcs 

In  this  section,  we  give  aji  efficient  algorithm  for  the  case  when  the  red  arcs  in  T 
are  non-intersecting  line  segments  which  form  the  boundary  of  a  simple  polygon  B 
(see  figure  1)  or,  more  generally,  non- intersecting  Jordan  arcs  whose  union  is  the 
boundary  of  a  simply-connected  region.  Let  A{r')  denote  the  arrangement  of  T' 
(i.e.  the  planax  map  induced  by  these  axes;  its  vertices  are  the  intersection  points 
of  the  axes  and  their  endpoints,  its  edges  axe  majcimal  connected  subsets  of  the  arcs 
which  do  not  contain  any  vertex,  and  its  faces  axe  mciximal  connected  subsets  of 
the  plane  which  do  not  contain  any  vertex  or  edge).  Let  z  denote  any  fixed  point 
on  dB.  The  following  Lemma  is  a  simple  but  a  key  observation  which  allows  us  to 
obtain  an  efficient  algorithm. 

Lemma  2.1  Let  T  be  the  face  of  A{T')  containing  the  point  z.  Then  dB  and  T' 
intersect  if  and  only  if  d!F  and  dB  intersect. 

Proof:  The  "if"  peirt  is  obvious.  For  the  "only  if"  paxt,  simply  follow  dB  from  z 
in,  say  clockwise  direction  until  the  first  intersection  with  T'  is  encountered,  which 
necessEirily  lies  on  d!F. 

The  above  lemma  suggests  that  we  do  not  have  to  look  at  the  entire  axrange- 
ment  of  T'  (whose  combinatorial  complexity  might  be  Q.{m^)  in  the  worst  case), 
but  instead  we  only  need  to  compute  a  single  face  of  A{T').  Guibas,  Shaxir  and 
Sifrony  [GSS88]  (see  also  [PSS88])  have  studied  the  complexity  of  a  single  face  in 
an  arrangement  of  Jordan  arcs.  They  have  shown 

Theorem  2.2  ([GSS88])  Given  a  collection  T'  =  {71,  72,  ...  ,7m}  of  Jordan  arcs 
in  the  plane  so  that  each  pair  of  them  intersect  in  at  most  s  points,  the  complexity 


of  a  single  face  of  A(T')  is  bounded  by  0(A,+2("^))  O'^d  it  can  be  computed  in 
0(Aa+2("^)log^  m)  time  (under  an  appropriate  m,odel  of  com.putation). 

Remark:  If  the  eircs  in  T'  axe  closed  (full)  Jordain  curves,  then  the  bound  on 
the  combinatorial  complexity  of  (resp.  the  time  to  compute)  a  single  face  of  their 
arrangement  can  be  replaced  by  0(\s{m))  (resp.  0(A,(m)log^  m))  (see  [SS87], 
[GSSS8]). 

The  combinatorial  bound  is  obtained  by  showing  that  the  circular  sequence  of 
arcs  in  the  order  in  which  they  appear  along  the  boundary  of  the  given  face  can 
be  written  as  a  Davenport-Schinzel  sequence  of  order  s  +  2.  The  algorithm  given 
by  [GSS88]  uses  a  divide  and  conquer  approach.  It  divides  T'  into  two  subsets 
r'l  =  {Yi,  ■■■  ,7m/2}  an«i  ^2  =  {l'm/2+i^  •••  '7m}  ^nd  recursively  computes  the 
face  ^i  of  ^(r'l)  and  the  face  T2  of  A{T'2)  containing  the  point  z.  The  desired 
face  ^  is  the  connected  component  of  ^i  n  ^2  which  contzdns  the  point  z.  Using  a 
relatively  simple  Hne-sweeping  algorithm,  a  subsequent  "red-blue  merge"  procedure 
then  obtains  ^  from  Ti,  T-i  in  0(Aj+2(m)logm)  time.  Hence,  the  overall  running 
time  is  0(A,+2("^)log'^  m). 

Remark:  Here  and  later  we  assume  a  model  of  computation  in  which  various 
basic  operations  involving  the  given  axes  (e.g.  finding  the  intersection  of  a  pair  of 
arcs,  or  of  an  axe  with  a  line,  or  testing  whether  a  point  lies  above  or  below  an 
axe,  etc.)  axe  assumed  to  require  0(1)  time.  We  eire  thus  more  interested  in  the 
combinatorial  complexity  than  in  the  algebraic  complexity  of  manipulating  such 
axes  (which  is  an  interesting  subject  for  further  study). 

Having  computed  T ^  all  we  have  to  do  is  test  whether  dT  and  dB  intersect. 
This  is  easily  done  by  a  variant  of  the  Bentley-Ottmzmn  algorithm  (see  [B079])  in 
time  0((A,+2(m)  +  n)  log(m  4-  n)).  Thus,  we  can  conclude 

Theorem  2.3  If  the  arcs  in  T  form,  the  boundary  of  a  simple  ("curved")  polygon 
B,  then  an  intersection  between  T  and  T'  can  be  detected  in  time 

O  (A,+2("^)  log^  m  +  (A,+2(m)  -|-  n)  log(m  +  n)) 
where  s  is  the  maximum  number  of  intersections  between  a  pair  of  arcs  in  T'. 
Proof:  Immediate  from  Theorem  2.2  and  the  discussion  given  above. 


Remarks: 

(i)  If  the  arcs  in  T'  are  also  line  segments,  then  we  can  use  the  ray-shooting 
technique  of  [CG85]  to  detect  an  intersection  of  a  blue  axe  with  the  red  polygon 
in  O(logn)  time,  after  0(n  log  log  n)  preprocessing. 

(ii)  If  the  axes  in  F'  axe  closed  Jordan  curves,  we  can  replace  5  -|-  2  by  5  in  the 
preceding  results. 

The  above  technique  can  be  extended  further  to  detect  an  intersection  between 
r  and  r'  when  each  of  A{T)  and  A{T')  is  a  connected  planar  graph.  Let  p  be  an 
endpoint  of  an  arc  in  T'  and  let  ^  be  the  face  of  >l(r)  containing  the  point  p.  Then 
we  have 

Lemma  2.4  F  and  F'  intersect  if  and  only  if  dT  and  F'  intersect  (see   figure  2). 

Proof:  The  "if"  paxt  is  obvious.  For  the  "only  if"  paxt,  let  5  be  a  red-blue  inter- 
section point.  Since  ^(F')  is  a  connected  graph,  there  exists  a  connected  path  IT 
from  p  to  q  along  the  edges  of  A{r').  Follow  II  from  p  until  its  first  intersection 
with  -4(F),  which  necessarily  lies  on  dT. 

Thus,  we  can  reduce  the  problem  of  detecting  an  intersection  between  F  and  F' 
to  the  problem  of  detecting  an  intersection  between  !F  (a  simply  connected  region) 
and  F',  which  we  can  solve  using  the  preceding  technique.  Therefore,  we  obtain 

Theorem  2.5  If  each  of  A{T),  A{T')  is  a  connected  planar  graph,  then  an  inter- 
section between  F  and  T'  can  be  detected  in  time  0(Aa+2(^T2-|-n)log^(77z-|-n)),  where 
s  is  the  maximum  number  of  intersections  between  a  pair  of  arcs  in  F  or  in  F'. 

Proof:  Using  Theorem  2.2,  we  can  compute  the  face  ^  in  time  0(A,+2(")log^  n). 
From  Lemma  2.4  it  follows  that  we  only  have  to  detect  an  intersection  between  d^ 
and  F',  which,  by  Theorem  2.3  and  the  following  remark,  can  be  done  in  time 

O  (A,+2(m)  log^  m  -I-  (A,+2(m)  +  n)  log(m  -|-  n)) 

and  hence  the  overall  rurming  time  is 

O  (\,+2im)  log^(m  +  n)  +  A,+2(n)  log^(m  -|-  n))       = 
O  (A,+2(m -I- n)log^(m -l-n))  . 

Remark:  Note  that  we  do  not  require  any  bound  on  the  number  of  intersections 
between  a  blue  axe  and  a  red  axe,  because  our  cdgorithms  stop  as  soon  as  one  such 
intersection  is  detected.  However,  the  algorithm  still  assumes  that  an  intersection 
between  a  blue  axe  and  a  red  arc  can  be  detected  in  0(1)  time. 


3      Applications  to  collision  detection  and  motion 
planning 

Using  the  intersection  detection  algorithm  of  the  previous  section,  we  present  an 
efficient  algorithm  for  the  following  problem: 

Given  a  simple  polygon  B  which  is  allowed  to  translate  and  rotate  in  the 
plane,  a  set  of  polygonal  obstacles  and  a  prescribed  continuous  motion 
ofB,  check  whether  B  collides  with  any  obstacle  during  this  motion. 

We  begin  by  reducing  this  problem  to  that  studied  in  the  previous  section. 
Let  O  =  {Oi,  O2,  ...  ,Oi}  denote  the  given  set  of  polygonal  obstacles  (having 
pairwise  disjoint  interiors).  We  can  represent  their  boundaries  as  a  collection  of 
"walls"  W  =  {Wi,  ...  ,  Wm}  where  each  Wi  is  an  edge  of  some  obstacle.  We  also 
define  two  coordinate  systems:  the  environment-frame  and  the  object-frame.  The 
environment-frame  is  the  usual  coordinate  system  in  which  the  obstacles  are  fixed 
and  B  moves.  On  the  other  hand,  the  object-frame  is  rigidly  attached  to  5,  so  in 
it  B  is  fixed  and  the  positions  of  the  obstacles  vary  depending  on  the  position  of  B 
(in  the  environment-frame).  The  position  of  B,  in  the  environment-frame,  can  be 
described  by  a  pair  Z  =  [.Y,  0],  where  X  =  (i,  y)  is  the  position  of  a  fixed  point  p 
in  B,  which  is  tadcen  to  be  the  origin  of  the  object-frame,  and  9  is  the  orientation  of 
the  X-axis  px  of  the  object-frame,  relative  to  the  environment-frame  (see  figure  3). 
Formally,  for  any  such  placement  Z  of  B  we  obtain  an  Euclidean  transformation  ^7 
which  maps  a  point  ^  in  the  object-frame  to  its  position  u  in  the  environment-frame 
as  follows. 

x{u!)     =     X -\- x{^)  ■  cosO  —  y(()  ■  sin9  (3.1) 

yiu)     =     y -{- x{0  ■  sine -{- yiO  ■  cose  (3.2) 

In  the  object-frame,  the  position  of  any  obstacle  O,  depends  on  the  position  Z  of  B 
in  the  environment-frame  and  is  obtained  by  the  corresponding  inverse  map  ^^  • 
Let  n  denote  a  continuous  motion  of  B,  given  as  a  continuous  path  11  :  [0, 1]  — > 

S"^  xS^.  We  call  11  a  basic  path  if  n([0, 1])  is  an  algebraic  curve  (in  x,  y  and  tan-) 

having  a  fixed  (and  small)  degree.  Examples  of  basic  paths  axe  translation  along 
a  straight  line  or  along  a  low-degree  algebraic  curve,  rotation,  sliding  with  a  fixed 
pair  of  vertices  of  B  touching  a  fixed  pair  of  walls,  etc.  (For  exposition  salce,  we 
keep  the  notion  of  a  bcisic  path  somewhat  loose  and  informal,  so  for  example  we  do 
not  specify  precisely  how  low  do  we  want  its  degree  to  be,  etc.)  We  assume  that 
n  is  either  a  basic  path  or  a  concatenation  of  basic  paths.   The  complexity  of  11, 
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denoted  by  |n|  =  A:,  is  the  number  of  basic  paths  from  which  it  is  composed,  and 
we  write  11  =  7ri||;r2||  •  •  •  \\TTk.  We  use  the  term  ^n  to  denote  the  map  which  sends 
each  point  ^  in  the  object-frame  to  its  orbit 

{*n(o(OUe[o,i]}. 

"i!^  is  defined  symmetrically  for  the  inverse  transformation  ^~'. 

Lemma  3.1  Let  q  be  a  fixed  point  in  the  environment-fraTne  and  let  B  move  along 
a  basic  path  U.  Then  ^nH?)  *■'  ^^^^  "-  basic  path  (albeit  its  maximum,  degree  m,ight 
be  higher).  Similarly,  if  W  is  a  wall,  then  '^^{W)  is  a  region  whose  boundary 
consists  of  a  constant  number  of  algebraic  curves  of  fixed  degree. 

Proof:  Let  q  =  (p,  a)  be  a  fixed  point  in  the  environment-frame  and  let  Tl{t)  = 
[(x(<),  y(i)),  B{t)]  denote  the  placement  of  B  at  time  i  as  it  moves  along  FI.  Suppose 
(^'  0  =  *n(t)(9),  then  by  (3.1)  and  (3.2) 


Therefore, 


p     =     x(t)  +  (,  ■  cosd{t)  -  C  •  sine{t) 
a     =     y{t)  +  (,  ■  sind{t)  +  C  •  cosdlt) 

^     =     (p  -  x{t))  ■  cos6{t)  +  (<7  -  y{t))  ■  sine{t) 
C     =     -{p  -  x{t))  ■  sine{t)  +  {a  -  y{t))  •  cosOit) 

Since  n  is  an  algebraic  curve  of  fixed  degree,  x{t),y{t),  sin6{t)  and  cos6{t)  (or 

61 1) 
equivalently,  tan——)  are  related  by  a  pair  of  polynomial  equations 


Q  ix,  y,  tan-j   =  0 

of  some  fixed  degree.  Eliminating  these  three  variables  from  the  four  equations  (cf. 
[VdW70]),  we  get  an  algebraic  curve  (in  ^,  Q  of  fixed  degree,  showing  that  ^n  (?) 
is  a  basic  path. 

Now  let  us  consider  the  area  swept  by  a  wall  W.  Without  loss  of  generality  let 
us  assume  that  W^  is  a  portion  of  the  y-axis.  It  is  easily  checked  that  any  point  on 
the  boundary  of  sE'n  (T^)  either  lies  on  one  of  the  paths  ^n  (9i)>  *n  (92),  where  qi, 
92  are  the  endpoints  of  H^,  or  lies  on  the  envelope  of  the  parametric  family  of  curves 
{*n  (^)  :  <  G  [0, 1]},  where  (  is  the  y-ax'is  (see  [Co36]  for  details).  For  each  t  e  [0, 1], 


let  n{t)  =  [x{t),y{t),e{t)].   Then,  the  points  {(,  Q  G  *n  (^)  (in  the  object  frame) 
satisfy  the  following  equations,  as  is  implied  by  equations  3.1,  3.2. 

/(^,  C,  0  =  3^(0  +  (cose{t)  -  ^sine{t)  =  0  (3.3) 

The  envelope  $  of  /(^,  (,  t)   =  0  (with  ^,  (  as  independent  variables  and  t  as 
a  parameter)  is  found  by  eliminating  t  from  the  pair  of  equations  /(^,  (,  t)  = 

— f{^i  Ci  0  —  0  (^^^  [Co36]).  In  our  case  these  equations  axe 
dt 

^cos6{t)  -  CsinOit)     =     -x{t) 

m 


^sind{t)  +  (cose{t)     = 


m 


(3.4) 


or 

x(t) 
^     =     -x{t)cos9{t)  +  -^sindit) 

x(t) 
C     =      x(t)sine{t)  +  -^^^cose{t) 

Q 

where  the  right-hand  sides  are  both  cdgebraic  in  x{t),  y{t),  tan-,  as  is  easily  checked. 
Eliminating  these  three  variables  as  was  done  in  the  case  of  a  point  shows  that  the 
envelope  $  is  algebraic  (of  degree  generally  larger  than  that  of  11,  but  nevertheless 
having  fixed  majcimum  value,  depending  on  the  degree  of  11).  Thus,  the  boundary 
of  ^^^(W^)  consists  of  a  constant  number  of  connected  subarcs  of  three  basic  paths 
^n  (91))  *^n  (52),  ^  (where  this  number  depends  on  the  pattern  of  intersections  of 
these  paths). 

D 

3.1      The  algorithm 

A  simple  approach  will  be  to  compute  the  region  swept  by  every  edge  e,  of  S  as 
B  moves  along  11  and  test  whether  it  intersects  any  wail  W.  Since  for  each  basic 
path  TTj,  the  region  swept  by  e^  is  of  fixed  complexity,  we  can  check  in  O^m)  time 
if  it  intersects  with  any  W.  There  are  k  basic  paths  and  n  edges,  so  the  total  time 
spent  will  be  0{k  ■  m  ■  n).  Using  the  aigorithm  of  Section  2,  we  show  that  we  can 
do  much  better  than  this  naive  approach. 

Without  loss  of  generality,  we  can  assume  that  the  initial  position  of  5  in  11  is 
a  free  position.  (In  practical  applications  this  can  be  expected  to  be  the  case;  in 
general,  we  can  run  a  standard  Hne  sweeping  algorithm  (such  as  that  of  [B079]), 
to  verify  that  B  starts  at  a  free  placement,  in  0{{n  +  m)log(n  -|-  m))  time.)   Let 
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7j  denote  the  boundary  of  the  region  TZ'j  in  the  object-fraxne,  swept  by  Wi  when  B 
moves  aJong  tt^,  i.e.  7^  =  9^~^(T^,),  emd  let  Ti  =  {7J,  75,  ...  ,7^}  (by  Lemma  3.1, 
each  7*  is  the  union  of  0(1)  algebraic  arcs). 

Lemma  3.2  B  intersects  a  wall  Wi  while  moving  along  TI  if  and  only  if  dB  inter- 
sects ^n  (W^i)  (in  the  object  frame).  Therefore  either  dB  intersects  an  arc  in  T,, 
or  dB  lies  inside    one  of  the  regions  TZy 

Proof:  We  only  prove  the  first  statement  of  the  Lemma  because  the  next  statement 
is  an  obvious  consequence.  Since  the  initial  position  of  B  is  free,  collision  between 
Wi  and  B  implies  that  there  exists  a  placement  Zo  6  11  such  that  dB  and  Wi 
intersect  at  some  point  C  €  "ifzoidB)  fl  Wi.  But  then  it  follows  by  definition  that 

Hence,  if  B  intersects  a  wall  W,  while  moving  along  IT,  then  dB  intersects  the  area 
swept  by  Wi  in  the  object  frame. 

a 

We  thus  first  test,  in  0{mk)  time,  whether  the  origin  p  lies  in  any  of  the  regions 
IZ'j.  If  so,  an  intersection  has  been  detected  and  we  can  stop  right  away.  If  not,  the 
previous  Lemma  implies  that  it  suffices  to  determine  whether  dB  intersects  an  arc 


in 


r  =  Ur.- 


1=1 


Theorem  3.3  Given  a  polygon  B  with  n  vertices,  a  path  U  =  7ri||7r2||  •  •  •  ||7rfc  and  a 
set  of  polygonal  obstacles  O  with  m  walls,  we  can  check  in  time 

O  (A,(mfc)log^(mfc)  +  (n  +  X,{mk))\og(n  +  mk)J 

if  B  collides  with  any  of  the  obstacles  in  O  when  it  moves  along  O.  Here  s  is  a 
constant  depending  on  the  maximum  algebraic  degree  of  the  basic  paths  of  II. 

Proof:  It  requires  only  0(mk)  time  to  check  if  the  origin  p  Hes  inside  any  of  the 
1Vj,  the  Eirea  swept  the  wall  Wi  when  B  moves  along  the  basic  path  n^.  We  can 
compute  r,  the  collection  of  the  boundciry  axes  of  all  the  regions  71'^  in  time  0{mk), 
and  Theorem  2.3  implies  that  we  can  detect  an  intersection  between  T  and  dB  in 
time 

O  (^X,(mk)  log^(mfc)  +  {n  +  X,{mk))  log(n  +  mkfj  , 

where  5  —  2  is  the  maximum  number  of  intersections  between  a  pair  of  arcs  in  T. 

a 


3.2      Applications  to  motion  planning 

In  this  section  we  describe  some  applications  of  the  above  coUision-detection  algo- 
rithm. This  algorithm  is  useful  when  mk  «  n,  because  in  this  case  it  takes  only 
0(A,(n)log^  n)  time  to  determine  if  the  polygon  B  colhdes  with  any  obstacle  in  O 
while  moving  along  the  prescribed  path  11.  Two  such  typical  cases  are:  (i)  There 
axe  0{n)  walls  in  the  environment  and  the  path  of  B  consists  of  0(1)  basic  paths, 
and  (ii)  The  environment  is  of  a  small  fixed  size,  i.e.  has  only  0(1)  walls,  but  the 
path  n  may  consist  of  as  many  as  0{n)  basic  paths. 

On  the  other  hand,  if  mfc  ^  n,  but  nk  ss  m,  and  /  =  0(1),  i.e.  there  axe  not  many 
obstacles  although  each  obstacle  can  have  many  edges,  then  we  can  use  a  similar 
approach  but  instead  of  working  in  the  object-frame,  we  use  the  environment -frame. 
The  basic  idea  is  as  follows.  Let  7?.^,  for  1  <  i  <  n,  1  <  _;'  <  /:,  denote  the  area 
swept  by  the  edge  e,  of  B,  when  B  moves  along  the  basic  path  tTj.  Let  7J  =  ^7^^  and 
Ti  =  {7J,  ...  ,71;}.    Using  Lemma  3.2  we  can  again  show  that  B  collides  with  aui 

n 

obstacle  Or  if  and  only  if  either  O^  lies  inside  any  of  TZ'^  or  dOr  intersects  T  =  [J  F,. 

1=1 
Let  rrir  be  the  number  of  edges  in  Or.  Since  it  takes  only  0{nk)  time  to  check  if  Or 

lies  inside  any  oiTt'j,  it  follows  from  Theorem  2.3  that  we  can  detect  an  intersection 

between  T  and  dOr  in  time 

O  \^\,+2{rik)  log^  nk  +  (m^  +  A,+2("^'))  log(n/:  +  mr))  , 

where  5  is  the  maximum  number  of  intersections  between  any  pair  of  zurcs  in  T. 
Summing  over  all  /  obstacles,  the  total  time  required  to  detect  a  collision  is 

O  (^\,+2{nk)\og^  nk  +  {m  +  \,+2{nk))\og{nk  +  mjj   =  0(A,+2(m)log^  m). 

because  nk  ^  m  and  /  =  0(1). 

The  mciin  advantage  of  either  of  these  two  vEiricmts  of  our  approach  is  therefore 
that  the  resulting  complexity  bound  does  not  involve  the  term  0{nin),  which  seems 
to  be  unavoidable  in  any  standard  motion  pleinning  algorithm  that  constructs  the 
explicit  configuration  space. 

Returning  to  our  apphcation,  there  axe  certain  motion  planning  problems  in 
which  one  CEin  define  a  single  canonical  motion,  and  show  that  if  there  exists  a 
collision-free  motion  from  the  initial  to  the  final  position,  then  this  canonical  path 
is  also  collision-free.  In  such  cases  we  can  apply  our  collision  detection  technique  to 
obtain  an  improved  motion  planning  algorithm.  Such  a  case  axises  in  Maddila  and 
Yap  [MY86],  who  have  studied  the  problem  of  moving  a  simple  n-sided  polygon  B 
axound  a  right-angle  corner  in  a  corridor.  They  have  presented  an  0{n'^)  algorithm 
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for  planning  such  a  motion;  using  our  techniques,  we  obtain  an  improved,  close  to 
linear  algorithm. 

A  corridor  with  one  right-angle  corner  is  an  infinite  L-shaped  region,  as  shown 
in  figure  3  (where  the  two  sides  of  the  corridor  extend  in  the  —x  direction  and  the 
4-y  direction,  respectively).  The  problem  of  moving  a  polygon  B  around  the  comer 
is  to  move  B  from  a  given  initial  position  Z/  =  [(i/,y/),  61]  in  the  left  side  of  the 
corridor,  to  a  final  position  Zp  =  [(iF,j/ir),  dp]  in  the  other  side  of  the  corridor, 
without  colliding  with  the  corridor  walls.  In  some  applications  Zp  may  not  be  given 
and  we  have  to  find  out  whether  B  can  be  moved  around  the  corner,  from  some 
given  initial  placement  Z/.  We  will  first  analyze  the  former  problem  and  then  show 
how  to  handle  the  latter  variant. 

It  is  shown  in  [MY86]  that  a  collision-free  motion  from  Z/  to  Zp  exists  if  and 
only  if  the  following  canonical  motion  of  B  is  collision-free: 

(i)  Translate  in  — y-direction,  until  B  touches  the  horizontal  wall  of  the  convex 
corner  of  the  corridor. 

(ii)  Translate  in  x-direction  until  B  touches  the  verticcd  wall  of  the  convex  corner. 

(iii)  Slide  B  while  touching  the  two  walls  of  the  convex  corner  until  the  orientation 
becomes  9p  (this  is  the  only  non-triviaJ  part  of  the  motion,  and  we  denote  it 

by  n). 

(iv)  Translate  in  y-direction  until  y{p)  =  yp. 

(v)  Translate  in  — i-direction  until  i(p)  =  xp. 

If  Z/  and  Zp  Eire  free  positions,  then  all  of  the  translational  motions  are  collision- 
free.  Thus  all  we  have  to  do  is  to  calculate  the  path  II  in  (iii)  and  verify  that  it 
is  collision-free.  This  is  done  by  applying  our  algorithm  (in  the  object  frame)  to 
detect  an  intersection  between  B  and  the  concave  corner  Eind  adjacent  walls  of  the 
corridor.  The  efficiency  of  the  algorithm  obviously  depends  on  the  complexity  of  II 
(i.e.  the  number  of  basic  paths  composing  it)  which  is  analyzed  as  follows. 

We  can  parameterize  II  using  the  orientation  of  B  as  it  changes  monotonically 
and  continuously  along  II  (see  [MY86]  for  details).  At  any  orientation,  the  rightmost 
vertex  u  oi  B  touches  the  verticzJ  wall  of  the  convex  corner  and  the  bottommost 
vertex  v  oiB  touches  the  horizontal  wall  of  the  corner.  Except  at  a  finite  number  of 
critical  orientations,  where  an  edge  adjacent  to  u  (resp.  v)  becomes  vertical  (resp. 
horizontal),  the  supporting  pair  {u,v)  is  unique  and  does  not  change  between  emy 
pair  of  consecutive  critical  orientations.  The  path  formed  by  sliding  B  with  a  fixed 
pair  of  vertices  touching  a  fixed  peiir  of  walls  is  known  as  a  glissette  [L06I];  each 
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point  on  B  traces  an  elliptic  arc,  and  the  path  (as  a  mapping  into  5^  x  5^)  is  easily 
seen  to  be  a  basic  path. 

Lemma  3.4  ([MY86])  The  sequence  (ui,  u,)  of  supporting  pairs  encountered  along 
n  changes  in  an  orderly  manner;  Ui,  u,  change  in  the  order  in  which  they  appear 
along  the  convex  hull  of  B.  Moreover,  the  number  of  supporting  pairs  is  at  most  2n 
and  they  can  be  found  in  0{n)  time. 

The  above  Lemma  implies  that  IT  consists  of  at  most  2n  basic  paths,  and  by 
Lemma  3.1  the  arcs  describing  the  "inverse"  motion  of  the  walls  in  the  object 
frame  of  B  are  algebraic  curves  of  some  fixed  low  degree.  Thus,  we  obtain 

Theorem  3.5  It  takes  0(A,(n)log^  n)  time,  (for  some  small  value  of  s  >  0)  to 
determine  whether  B  can  be  moved  around  a  right-angle  corner  in  a  corridor. 

Proof:  Immediate  from  Lemma  3.4  and  Theorem  3.3,  noting  that  m  is  constant 

and  k  <  2n. 

D 

If  Zp  is  not  specified,  then  moving  B  around  the  corner  amounts  to  sliding  it  as 
in  the  above  canonical  motion  until  it  reaches  an  orientation  Op  at  which  it  can  be 
translated  vertically  into  the  upper  portion  of  the  corridor,  for  which  it  is  necessary 
and  sufficient  that  the  width  of  B  in  the  direction  perpendicular  to  dp  should  be 
less  than  the  width  of  the  upper  portion  of  the  corridor.  This  observation  enables 
us  to  find  the  smallest  such  9p  in  0{n)  time  using  the  standard  "rotating  calipers" 
method.  Having  found  6p,  the  problem  reduces  to  the  one  just  studied,  and  the 
same  algorithm  can  be  applied. 

We  can  extend  our  algorithm  to  a  more  general  corridor  in  which  the  exterior 
angle  of  the  concave  corner  is  not  less  than  the  interior  angle  of  the  convex  corner 
(see  figure  4a).  Indeed,  our  algorithm  relies  on  the  ability  to  reduce  any  collision- 
free  motion  of  B  around  the  corner  to  a  canonical  one  in  which  B  slides  along  the 
walls  of  the  lower  (convex)  portion  of  the  corridor.  The  conditions  eissumed  axe 
easily  seen  to  imply  that  a  "retraction"  similar  to  that  used  above  (first  translate 
B  to  the  right  until  it  touches  the  right  wall,  and  then  slide  it  downwards  along  the 
wall  until  the  lowest  vertex  of  B  touches  the  lower  horizontal  wall  of  the  corridor) 
can  be  applied  in  this  case  (see  figure  4b  for  an  example  where  such  a  retraction 
is  not  continuous  when  the  above  conditions  axe  not  satisfied).  Once  again,  we  can 
show  that  this  canonical  motion  consists  of  at  most  2n  basic  paths,  leading  to  an 
algorithm  similar  to  that  presented  above. 

If  S  is  a  star-shaped  polygon,  we  can  somewhat  improve  the  algorithm,  because 
we  can  replace  the  calculation  of  a  single  face  in  A{T)  (as  in  [GSS88])  by  the  calcu- 
lation of  a  certain  lower  envelope  of  F  (where  T  is  the  collection  of  arcs  describing 
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the  inverse  motions  of  the  walls  as  above).  To  this  end,  let  ^  be  a  point  in  the  kernel 
of  5,  and  consider  each  axe  7  of  F  as  a  function  r  =  ^{6)  in  polar  coordinates  about 
^  in  the  object  frame  (if  necessary,  we  break  such  a  7  into  a  constant  number  of 
pieces,  each  monotone  with  respect  to  9).  The  lower  envelope  Mr  of  T  in  these 
coordinate  is  defined  as 

MriO)  =  min{7(^)} 

Lemma  3.6  B  collides  with  a  wall  forming  the  concave  corner  of  the  corridor  while 
moving  along  the  canonical  path  if  and  only  if  B  intersects  the  lower  envelope  of 
A^r  defined  above. 

Proof:  Lemma  3.2  imphes  that  B  collides  with  a  wall  VT,  while  moving  along  the 
canonical  path  11  if  and  only  if  B  and  ^^{Wi)  intersect.  Therefore,  intersection  of 
Mt  and  B  implies  that  B  collides  with  some  Wi  while  moving  along  11. 

Conversely,  let  ^  G  5  n  d'^^{Wi).  ^  is  a  point  lying  on  one  of  the  curves  in  F, 
and  if  ^  ^  Mt-,  then  some  other  point  a  on  the  segment  C,^  must  he  on  this  lower 
envelope.  Since  B  is  stax-shaped  and  C,  lies  in  its  kernel,  the  entire  segment  (,  (,  and 
thus  also  <T  lies  inside  B.  Therefore,  Mr  and  B  intersect. 

An  intersection  between  My  and  B  can  be  easily  detected  as  follows.  Sort  the 
"breakpoints"  of  A^r  in  anguleir  order  about  ^  (where  each  breakpoint  is  either  an 
endpoint  of  some  arc  7  €  F,  a  point  of  radial  tangency  on  such  an  axe,  or  a  point 
where  A^p  is  simultaneously  attained  by  two  arcs  in  F),  and  similarly  sort  the 
vertices  of  B  in  angulax  order  about  (^  (note  that  the  latter  sorted  list  is  the  same  as 
the  order  of  the  vertices  of  B  along  its  boundary).  Merge  these  two  lists  to  obtain 
an  overall  list  C  of  orientations  about  (^,  with  the  property  that  for  each  interval 
2  between  any  pair  of  consecutive  orientations  in  £,  Mt  is  attained  over  Z  by  a 
single  Eire  7j  cind  the  portion  of  dB  seen  from  C  in  directions  in  2"  is  (a  portion  of)  a 
single  edge  ej.  We  can  thus  detect  in  constant  time  an  intersection  between  B  and 
Mt  within  the  angular  sector  defined  by  J.  There  axe  0(Aj(n))  breaJcpoints  in  A^r 
(for  some  small  fixed  5  >  0),  so  the  length  of  C  and  the  time  required  to  detect  an 
intersection  between  Mt  and  B,  axe  both  0{\3{n)).  Mt  can  be  calculated  (in  the 
required  sorted  order)  in  time  0(A,(n)logn),  using  a  standard  divide  and  conquer 
technique  [At85].  We  thus  have 

Theorem  3.7  It  takes  only  0{X,{n)\ogn)  time  (for  some  small  s  >  0)  to  deter- 
mine whether  a  star-shaped  polygon  B  can  be  moved  around  a  corner  in  a  corridor. 

Remark:  A  similar  technique  applies  if  S  is  a  monotone  polygon,  in  which  case 
the  radial  lower  envelope  is  replaced  by  appropriate  lower  and  upper  envelopes  in 
the  direction  of  monotonicity  oi  B;  we  leave  details  of  this  easy  extension  to  the 
reader. 
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4     Red  blue  intersection  detection  in  general 

If  the  arcs  in  F  or  in  T'  do  not  form  a  simply  connected  region  (or  worse  still  they  axe 
arbitrarily  intersecting  and  their  axraxigements  are  disconnected),  then  Lemma  2.1 
does  not  hold  and  we  may  have  to  search  for  an  intersection  in  more  than  one  face 
of  A{T').  Two  difficulties  can  then  arise.  One  is  that  the  complexity  of  many  faces 
in  such  an  arrangement  can  be  much  higher  than  lineeir.  Second,  the  algorithm  of 
[GSS88]  does  not  generalize  to  many  faces,  at  least  not  for  general  axes. 

As  discussed  earlier,  a  sweep-line  approach  is  doomed  to  be  quadratic  in  the 
worst  case,  because  it  may  have  to  sweep  across  many  red-red  or  blue-blue  inter- 
sections. A  technique  of  Chazelle  for  point  location  in  algebrciic  manifolds  [Ch85] 
can  be  used  (in  the  case  of  sdgebraic  arcs)  to  obtain  a  slightly  subquadratic  algo- 
rithm, but  otherwise  we  are  not  aware  of  einy  previous  algorithm  which  achieves 
(substantially)  subquadratic  performance  for  this  general  red-blue  intersection  de- 
tection problem.  In  this  section  we  present  algorithms  of  this  kind.  We  first  describe 
a  deterministic  algorithm  which  works  for  general  axes,  and  then  obtain  in  several 
restricted  but  useful  special  cases  more  efficient  randomized  algorithms. 

4.1      Deterministic  algorithm 

Our  deterministic  algorithm  consists  of  three  phases.  In  the  first  phase,  we  paxtition 

TTl 

T'  into  £i  <  —  sets  F'j,  Fj,  ...  ,F/  ,  for  some  parameteT  rj  to  be  chosen  later,  so 

that  each  F^,  for  1  <  j  <  ii,  has  rrij  arcs,  where  rj  <  m^  <  2ri  and  A{T'j)  forms 
a  connected  planax  graph.  F^^  (which  may  be  empty)  contains  the  remaining  arcs 
of  F',  so  that  each  connected  component  of  A{T(^)  (i.e.  of  the  union  of  the  axes 
in  T'f^  ^  )  has  less  than  rj  aires.  For  each  1  <  j  <  ^i,  we  check  whether  an  axe 
of  T'j  intersects  an  axe  of  F,  using  a  modified  version  of  the  technique  of  Section 

Tt 

2  (to  be  described  below).    In  the  second  phase  we  partition  F  into  ^2  <  —  sets, 

Fi,  F2,  . . .  ,  F/j ,  satisfying  conditions  analogous  to  those  of  the  partition  of  F',  and 
for  each  1  <  j  <  ^2,  cheek  whether  Fj  intersects  F'.  Finally,  in  the  third  phase 
we  check  whether  F^^  and  F^^  intersect.  Together,  the  three  phases  will  detect  an 
intersection  between  F  and  F',  if  one  exists. 

The  first  phase  of  the  algorithm  proceeds  as  follows.  We  paxtition  F'  into 
Fj,  ...  ,FJj  using  a  line  sweep  approach.  (As  above,  we  assume  all  arcs  to  be 
x-monotone.)  We  sweep  F'  from  left  to  right  by  a  vertical  line;  every  time  we  sweep 

^In  what  follows  we  use  the  term  "connected  component"  of  an  arrangement  to  refer  to  a  compo- 
nent of  the  union  of  its  arcs,  whereas  the  term  "face"  will  continue  to  refer  to  a  connected  component 
of  the  complement  of  the  union  of  these  arcs. 
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through  an  intersection  of  two  arcs,  we  merge  the  connected  components  of  A{T') 
containing  them,  eind  whenever  a  component  is  found  to  contain  more  than  ri  arcs, 
we  remove  all  these  arcs  from  F'.  In  more  detail,  as  in  any  standard  line  sweep  al- 
gorithm, we  maintain  a  Hst,  denoted  by  Q,  of  arcs  currently  intersecting  the  sweep 
line,  sorted  in  increasing  order  of  y  coordinate  of  their  intersections.  In  addition 
we  maintain  the  connected  components  of  A{T'  fl  h~),  where  T'  is  the  set  of  axes 
that  have  been  encountered  by  the  sweep  line  so  far  but  have  not  yet  been  deleted 
by  the  edgorithm,  and  h~  is  the  half  plane  lying  to  the  left  of  the  sweep  line.  For 
each  axe  7'  6  F',  we  store  the  connected  component  ^(7')  which  contains  it  (if  an 
arc  7'  has  not  been  encountered  by  the  sweep  line,  0(7')  is  not  defined).  At  every 
event  point,  during  the  sweep  (i.e.  an  endpoint  of  one  arc  or  an  intersection  point 
between  two  eircs  in  F'),  apart  from  the  stEmdeird  list  and  priority  queue  updating 
operations,  we  perform  the  following  additional  tasks. 

(i)  If  we  reach  a  left  endpoint  of  some  arc  7',  then  we  create  a  new  connected 
component  C{Y)  containing  only  7'  (sweeping  through  the  right  endpoint  of 
an  axe  requires  no  special  action). 

(ii)  If  we  reach  an  intersection  point  of  two  aires  7',  7"  and  C{f')  is  different 
from  C(j"),  then  we  merge  the  two  connected  components  containing  7',  7" 
respectively. 

(iii)  If  suiy  connected  component  C  is  found  to  have  more  than  rj  eires,  then  we 
output  it  as  a  set  F^,  and  delete  C  and  all  of  its  axes  from  F'  (and  from  Q). 
(Note  that  C  need  not  be  a  full  component  of  -4(F'),  but  may  be  only  a  subset 
of  such  a  component.)  We  also  remove  all  endpoints  and  intersection  points 
along  these  axes  from  the  priority  queue.  The  priority  queue  is  also  updated 
by  inserting  into  it  all  the  intersection  points  between  newly  adjacent  pairs  of 
arcs  in  Q. 

It  is  easy  to  see  that  the  number  of  axes  in  each  connected  component  processed 
during  the  sweep  is  between  rj  and  2ri .  The  axes  not  deleted  by  the  sweep  are  put 
into  the  set  FJ^ . 

Next,  for  each  I  <  j  <  £\,  v/e  check  whether  F^  intersects  F  using  the  following 
procedure. 

We  first  construct  the  axreingement  ^(F^)  of  the  arcs  in  F^,  using  a  line  sweep 
technique.  If  the  two  endpoints  of  an  arc  7  €  F  he  in  two  different  faces  of  A{T'j), 
then  obviously  7  intersects  some  axe  of  F^.  If  they  lie  in  the  same  face  fi  of  A{T'j), 
then  7  need  not  intersect  F^;  however  if  it  does,  it  has  to  intersect  a  subaxe  lying 
on  the  boundary  of  /,  (as  in  Lemma  2.1).  Therefore,  the  next  step  is  to  determine 
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for  each  axe  7  G  F,  which  faces  contain  the  endpoints  of  7.  This  can  be  easily 
done  during  the  hne  sweeping  procedure  which  produces  A{T'^)  (see  also  [Pr79]),  in 
overall  time  0{{rl  +  n)log(ri  +  n)).  To  reiterate,  our  sweep  performs  the  following 
steps  for  each  7  €  F. 

(i)  If  the  two  endpoints  of  7  lie  in  two  different  faces  of  ^(F^),  then  we  have 
found  a  red-blue  intersection  and  we  stop  right  away. 

(ii)  If  the  two  endpoints  of  7  lie  in  the  same  face  /,  of  A{r'j),  then  we  assign  7  to 

Thus,  if  the  sweep  does  not  detect  an  intersection,  it  produces  a  partition  F/, ,  . . .  ,  F/^ 
of  F,  where  p  =  0{rl),  such  that  all  eircs  in  F/,  have  both  of  their  endpoints  in  the 
face  /,  of  the  arrangement.  Let  n,  =  |F/,  |.  Since  every  arc  of  F  is  in  exactly  one 
F/, ,  IZf=i  n.  =  n.  As  argued  above,  it  suffices  to  check  the  arcs  of  each  F/,  for 
intersection  with  the  edges  of  /,.  A  simple  way  of  doing  it  is  to  take  an  arc  of  F/, 
and  check  its  intersection  with  all  edges  of  /,.  But  in  the  worst  case  /,  may  be 
bounded  by  f2(ri )  arcs  amd  all  arcs  of  F  may  lie  in  this  face,  in  which  case  this  naive 
procedure  will  be  too  expensive.  The  following  Lemma  suggests  a  way  to  improve 
this  procedure  by  exploiting  the  property  that  the  boundary  of  each  face  in  A{T'j) 
is  connected. 

Lemma  4.1  Let  fi  be  a  (simply  connected)  face  of  A{T'j),  and  let  T^  be  the  un- 
bounded face  of  A{Tf^).  An  arc  7  6  F/,  intersects  F^  if  and  only  if  dTi  and  df, 
intersect  (see  figure  5).  (For  unbounded  faces  fi  of  A{Tj),  take  Tx  to  be  the  face  of 
A{Tj,)  that  contains  some  point  of  df;.) 

Proof:  Using  the  argument  of  Lemma  2.1  we  can  prove  that  F/,  and  F^  intersect 
if  and  only  if  5/,  and  F/,  intersect.  By  assumption,  5/,  is  a  simple  closed  Jordan 
curve,  or  a  simple  unbounded  Jordan  axe.  Talce  any  point  2  6  5/,,  and  let  Ti  be 
the  face  of  A{T j^)  that  contain  z  (it  is  easily  verified  that  if  /,  is  bounded  then  T, 
is  the  unbounded  face  of  >t(F/,)).  The  claim  now  follows  from  another  application 

of  Lemma  2.1,  in  which  the  roles  of  F  and  F'  axe  interchanged. 

D 

We  thus  compute  the  unbounded  faces  !Fi  for  each  F/,  and  then  can  easily  detect 
an  intersection  between  any  3/,,  d^i  by  performing  another  line  sweep. 

This  completes  the  description  of  the  first  phase  of  our  algorithm.  The  second 
phase  proceeds  in  a  completely  symmetric  manner  (where  we  paxtition  F  into  con- 
nected components  of  size  r-^  (to  be  specified  later)  instead  of  ri).  For  the  final 
phase,  we  first  construct  A{T'(^)  and  A{Ti^).  Then  by  sweeping  the  two  arrange- 
ments with  a  verticed  hne  we  check  whether  an  edge  of  A{T[^ )  intersects  an  edge 
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of  -4(r<2 )  (in  this  finaJ  sweep  we  are  willing  to  encounter  all  red-red  and  blue-blue 
intersections  as  well). 

The  correctness  of  the  algorithm  follows  easily  from  the  analysis  given  above. 
As  to  Einalysis  of  its  time  complexity,  we  first  need  the  following  simple  lemma. 

Lemma  4.2  //  each  connected  component  of  A{T)  has  at  most  r  arcs,  then  A{T) 
has  only  0{nr)  vertices. 

Proof:  Let  Fi,  ...  ,Tk  denote  the  sets  of  8ircs  constituting  each  connected  compo- 

r 
nent  of  ^(F).  If  two  sets  F,,  Fj  have  at  most  -  arcs  each,  we  merge  them  together. 

We  repeat  this  process  until  there  are  no  two  such  subsets.  Then  each  F,  (except 

r 
possibly  one)  has  at  least  -  but  less  than  r  arcs,  so  the  number  of  these  sets  are  at 

2n 

most hi  and  each  A{Ti)  has  0{r'^)  vertices.  Since  A{Ti)  and  A{Tj),  for  j  /  i, 

r 

h  1 )  X  0{r^)  =  0{nr)  vertices  in  >t(F). 

Lemma  4.3  Partitioning  F'  into  T\,  ...  ,F^  in  the  first  phase  of  the  algorithm 
can  be  accomplished  in  O(mrilogm)  time. 

Proof:  Inspecting  the  line  sweeping  cdgorithm  used  to  produce  the  partitioning, 
we  see  that  most  of  its  steps  can  be  performed  quite  efficiently.  For  example,  main- 
taining the  connected  components  of  .4(F')  (creating  new  connected  components 
and  merging  pairs  of  components)  can  be  easily  done  in  time  0{ma{m)),  using  the 
union-find  data  structure  [Ta83].  Also,  the  time  needed  to  delete  arcs  from  Q  every 
time  a  component  is  produced  is  0{m  log  m)  because  each  arc  is  deleted  only  once. 
The  most  expensive  part  of  the  eJgorithm  is  the  maintenance  of  the  priority  queue, 
which  teikes  time  0{K  logm),  where  K  is  the  number  of  endpoints  and  intersection 
points  encountered  by  the  algorithm.  Therefore,  to  prove  the  Lemma,  we  only  have 
to  show  that  K  =  O(mri). 

Let  r[,  . . .  ,T'f  be  the  connected  components  produced  by  the  algorithm.  Since 
the  algorithm  merges  connected  components  whenever  it  sweeps  through  an  inter- 
section point  of  two  Eircs  lying  in  two  different  components,  it  never  encounters  an 
intersection  point  between  axes  of  two  different  F'^  and  F^  (more  precisely,  it  may 
have  added  such  points  into  the  priority  queue,  but  it  never  reaches  these  points 
during  the  sweep).  Moreover,  at  every  intersection  swept  by  the  line,  only  a  con- 
stant number  of  events  are  added  to  the  priority  queue.  Hence,  the  total  number 
K  of  intersections  encountered  by  the  algorithm  is  proportional  to  the  number  of 
intersections  swept  through,  which  in  turn  is  at  most  the  total  number  of  vertices 
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in  all  the  subaxrangements  A{T'j),  I  <  j  <  ii-  For  j  <  £i,  T'^  has  less  than  2ri  arcs, 
therefore  A{T'.)  has  0(rj)  vertices.  Furthermore,  since  each  connected  component 
of  A{T'f  )  hsis  less  thaji  ri  arcs,  Lemma  4.2  implies  that  this  subaxrangement  has 

only  O(mri)  vertices.  Since  ^i  <    —  ,  the  total  number  of  intersection  points  K  is 

O(mri).  Hence,  it  takes  only  0(mri  logm)  time  to  partition  T'  into  r[,  ...  ,r^, . 

D 


Lemma  4.4   The  time  spent  in  the  first  phase  of  the  algorithm  is  bounded  by 

0{mri  log(ri  +  n)  +  — A,+2(n)  log^  n). 

Proof:  It  follows  from  the  previous  Lemma  that  T[,  ...  ,T'(^  can  be  obtained  in 
time  0{mri  logm).  Therefore,  we  only  have  to  bound  the  time  spent  in  detecting 
an  intersection  between  F^,  for  1  <  j  <  ^i  and  F. 

If  |F/, I  =  n,,  then  Lemma  2.2  imphes  that  d^,  can  be  computed  in  time 
0(A,+2("«)log^  ":)•  By  using  a  line  sweep  technique  we  can  determine  in  time 
0{{mi  +  Aa+2(n,))log(n,  +  m,))  whether  dTi  £md  fi  intersect,  where  m,  is  the  num- 
ber of  edges  in  dfi.  Summing  these  costs  over  all  faces  fi  of  A{r'j),  we  get 

p 

j^O(m,log(n, +  mi)  +  A,+2(".)log^  n.) 
1=1 

=     Oi^m,  log(n,  +  m.)  +  J^  A,+2(n.) log^  ".• ) 
\i=i  i=\  ) 

The  total  number  of  edges  in  ^(F^)  is  bounded  by  0{t\),  therefore  IIf=i  m,  = 
0{r\).  Since  XZf=i  "i  =  ",  and  ^(F^)  cmd  F/j,  ...  ,F/p  can  be  computed  in  time 
0{t\  logTi)  eind  0((r^ +  n)  log(ri +n))  respectively,  the  total  time  spent  in  detecting 
an  intersection  between  F  and  F^  is 


Since,  t\  < 


m 


O  (r^  log(ri  +  n)  +  A3+2(n)log^  ") 


,  the  overall  time  spent  in  the  first  phase  is  at  most 


O  {—r\  ■  log(n  +  ri )  +  -A,+2(n)  •  log^  "  ) 
=     O  (mri  log(n  +  ri)  +  —  •  A,+2(n)log^  nj 
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Using  a  symmetric  analysis,  it  follows  that  the  second  phase  requires 
O  (n  ■  r2log(m  +  r2)  -\ A,+2(m)log^  m  j 

time.  As  to  the  third  phase,  Lemma  4.2  implies  that  A{T'f^)  has  only  0{mri) 
vertices,  therefore  we  can  easily  compute  it  in  time  0{mri  logm).  Similarly,  A{Ti^) 
has  0(nr2)  vertices  aind  cam  be  computed  in  time  0(nr2  log  n).  Once  we  have  A{T'(  ) 
and  A{Tt^),  we  can  easily  check  in  time  0((mri  +  nr2)log(m  +  n))  whether  they 
intersect,  which  bounds  the  running  time  of  the  third  phase.  We  thus  obtain  the 
main  result  of  this  section 

Theorem  4.5  Given  a  set  of  n  "red"  Jordan  arcs  T,  and  another  set  of  m  "blue" 
Jordan  arcs  T' ,  we  can  detect  in  time 

O  ((niyj \,+2in)  +  n^J \,+2im))\og^-^{m  +  n) 

whether  T  and  T'  intersect,  where  s  is  the  maximum  number  of  intersections  between 
a  pair  of  arcs  in  T  or  in  T'. 

Proof:  Summing  the  cost  of  all  the  three  phases  of  the  algorithm  and  choosing 
''i  =  y^3+2{''^)^og{m  +  n),  r2  =  \/Aj+2(m)log(m  +  n),  the  overall  running  time 
becomes 

O  ((myJX,+2{n)  +  n^A,+2(m))log^-^(m  +  n)) 

Remarks:  In  some  special  cases  we  can  do  somewhat  better,  particularly  when 
n  '^  m.  We  mention  two  such  cases. 

(i)  If  r'  can  be  partitioned  into  subsets  of  the  appropriate  size  r  so  that  the  bound- 
ary of  each  face  in  the  arrangement  of  each  subset  is  simply  connected,  then  we 
do  not  need  the  last  two  phases  of  the  above  algorithm.  Therefore,  it  easily  fol- 
lows that  the  running  time  of  the  algorithm  becomes  0{mJXa+2{n)^og^'^  n), 
which  is  certainly  better  than  the  running  time  of  the  above  algorithm  for 
n  ':^  m.  A  case  in  which  we  can  ensure  that  each  /,  is  simply  connected  is 
when  the  arcs  in  T'  axe  concatenated  to  one  another  at  their  endpoints  to  form 
a  single  connected  (possibly  self-intersecting)  chain.  This  situation  typically 
arises  in  collision-detection  problems  of  the  sort  studied  in  Section  2  and  3. 

(ii)  Another  special  case  in  which  we  can  get  a  slightly  improved  bound  is  when 
the  axes  in  T  are  non-intersecting.    Again,  we  do  not  need  last  two  phases. 
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We  partition  T'  arbitrarily  into     —     sets  T',  each  of  size  at  most  ri,  and 

then,  for  each  T^,  as  in  the  first  phase  of  the  general  algorithm,  we  compute 
F/, ,  ...  ,  F/  .  Since  arcs  in  F  are  non-intersecting  we  do  not  have  to  compute 
any  specific  face  of  F/, ,  but  instead  we  can  detect  an  intersection  between  5/, 
and  F/,  by  performing  a  line  sweep.  By  choosing  ri  =  y/n,  one  can  easily 
check  that  the  overall  running  time  is  bounded  by  0{my/n\ogn). 

(iii)  Again,  if  all  the  arcs  in  F,  F'  are  closed  Jordan  curves,  then  we  can  replace 
s  -t-  2  by  5  in  the  preceding  results. 

4.2      Reporting  red-blue  intersections 

If  all  the  arcs  in  F,  F'  are  hne  segments,  then  we  can  extend  our  algorithm  to 
report  all  red-blue  intersections  in  time  that  depends  on  the  actual  number  of  such 
intersections.  Again,  the  algorithm  consists  of  three  phases.  We  partition  F  into 
subsets  Fj,  for  1  <  j  <  £2,  and  F'  into  subsets  F^,  for  1  <  j  <  ^1,  as  described 
above.  In  the  first  phase  we  report  the  intersections  between  F  ajid  (F'  —  F^^),  in 
the  second  phase  we  report  the  intersections  between  F^^  and  (F  —  F/^),  and  finsdly 
in  the  third  phase  we  report  the  intersections  between  F^^  and  F^j . 

The  first  pheise  proceeds  as  follows.  As  in  the  detection  algorithm,  we  consider 
each  F^  separately.  Each  face  of  A{r'j)  is  a  simple  polygon  (or  "semi-simple"  in  the 
sense  that  its  boundary  may  touch  itself,  when  two  sides  of  a  segment  bound  the 
same  face,  but  not  cross  itself).  We  preprocess  each  face  of  A{r'j)  for  answering  ray- 
shooting  queries  using  the  techniques  of  [CG85],  [GHLST87]  (which  also  apply  to 
semi-simple  polygons).  Let  p,  a  be  the  endpoints  of  a  segment  7  G  F,  and  let  /p,  /„ 
be  the  faces  of  A{T'j)  containing  these  endpoints  (as  determined  in  the  first  phase  of 
the  detection  algorithm).  By  performing  a  ray-shooting  query  we  find  out  the  edge 
e  of  fp  which  pa  hits  first.  If  there  is  no  such  edge,  then  7  does  not  have  any  red-blue 
intersection;  otherwise  the  intersection  point  t]  is  obviously  a  red-blue  intersection 
and  we  report  it.  Let  /^  be  the  other  face  of  ^(F^)  adjacent  to  e.  We  now  repeat 
the  same  process  for  the  segment  rfa  in  /,,.  We  continue  this  until  we  reach  the 
other  endpoint  a  of  7.  Since  for  each  axe  7  G  F  all  ray-shooting  queries  (except 
the  last  one)  return  distinct  red-blue  intersections,  we  can  charge  the  Icist  query 
to  7  and  the  remEiining  queries  to  the  intersections.  The  preprocessing  of  face  /, 
requires  0(|/,|  log  \fi\)  time  (or  0(|/,|loglog  |/, |)  using  a  more  involved  algorithm  of 
[GHLST87])  and  a  ray-shooting  query  cjin  be  einswered  in  time  O(logri),  therefore 
all  k'^  red-blue  intersections  between  F  and  F'  can  be  reported  in  time  0(n  log  n  + 
Tj  logri  -|-  fcj  logri).  Using  the  same  argument  as  in  Lemma  4.4  we  can  show  that 
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the  total  time  spent  in  the  first  pheise  is  0(mri  log(ri  +  n)  +  nlogn  +  fcilogri), 
where  ki  is  the  number  of  red-blue  intersections  between  F  and  F'  —  F^  .  Similaxly, 
all  k2  red-blue  intersections  between  F^^  and  F  —  T(^  can  be  reported  in  the  second 
phase  in  time  0(nr2  log(r2  +  m)  +  m  log  m  +  k2  log  r2 ).  Since  in  the  third  phase  we 
can  afford  the  calculation  of  red-red  and  blue-blue  intersections,  we  Ccin  report  all 
red-blue  intersections  by  sweeping  over  all  the  arcs  in  T'(  ,  Te^.  Thus,  summing  up 
all  the  costs  and  choosing  ri  =  y/n  and  r2  =  y/rn,  we  obtain 

Theorem  4.6  //  all  the  arcs  in  T,  T'  are  line  segments,  then  all  k  red-blue  inter- 
sections can  be  reported  in  time  0({my/n  +  riy/m  +  A:)log(Tn  -|-  n)). 

Remark:  The  above  algorithm  also  shows  that  if  all  arcs  are  line  segments,  then 
the  red-blue  intersection  detection  problem  can  be  solved  in  deterministic  time 
0{{my/n  +  riy/rn)  log(m  -f  n)),  which  is  slightly  better  than  the  running  time  of  the 
general  algorithm. 

4.3      Relation  between  many  faces  and  red-blue  intersection 
detection 

In  this  subsection  we  explore  the  close  relationship  between  the  problem  of  comput- 
ing many  faces  in  the  arrangement  of  a  given  set  of  Jordem  axes  and  the  red-blue 
intersection  detection  problem  for  those  arcs.  That  is,  we  show  that  an  efficient 
edgorithm  for  the  former  problem  can  be  used  to  obteiin  an  efficient  algorithm  for 
the  latter  one.  The  basic  idea  is  that  a  red-blue  intersection  can  be  detected  by 
knowing  only  a  few  faces  in  A{T)  and  A{T'),  which  we  prove  in  the  following  lemma. 

Lemma  4.7  Let  T  =  {dfi,df2,  . . .  ,dfk}  (resp.  T  =  {df[,df'^,  . . .  ,df'k.])  be  the 
collection  of  the  (connected  components  of  the)  boundaries  of  the  faces  in  A{T) 
(resp.  A{T'))  which  are  either  unbounded  or  contain  at  least  one  endpoint  of  (a 
bounded)  arc  either  in  F  or  V .  Then  F  and  F'  intersect  if  and  only  if  T  and  T' 
intersect. 

Proof:  The  "if"  part  is  obvious.  For  the  "only  if"  part,  suppose  first  that  all  arcs  in 
r,  F'  are  bounded.  Let  a  be  a  red-blue  intersection  point.  Let  >^(F^)  (resp.  .4(F^)) 
denote  the  connected  component  of  >t(F)  (resp.  >l(F'))  containing  the  point  a  (see 
figure  6).  It  is  obvious  that  the  rightmost  point  r'  of  >1(F^)  lies  in  its  unbounded 
face.  Therefore  the  boundary  df[,  of  the  unbounded  face  of  ^(F^)  belongs  to  the 
set  T' .  Let  r  be  the  rightmost  endpoint  of  >l(Fj,),  which  again  must  lie  in  the 
unbounded  face  fr  of  Ai^^).,  so  that  dfr  is  in  T .  Without  loss  of  generality  assume 
that  r  is  to  the  right  of  r',  which  implies  that  r  lies  in  f[,.  There  exists  a  (connected) 
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path  TT  in  A{r„)  from  r  to  a,  and  let  a*  be  the  first  red-blue  intersection  between 
A{Tff)  and  A{T'^)  along  n  (note  that  tt  could  intersect  other  components  of  A{T') 
as  well,  but  we  ignore  such  intersections).  Clearly,  a*  must  lie  on  df^,.  If  a*  also 
lies  on  dfr ,  then  dfr  and  3/^,  intersect  which  implies  the  assertion  of  the  lemma. 
Otherwise  we  follow  dfl,  from  r'  to  cr*,  and  let  <7**  be  the  first  red-blue  intersection 
along  this  boundary.  Since  a**  hes  on  df^,,  and  also  on  the  boundary  of  the  face  of 
A{T)  containing  r',  we  again  obtain  the  assertion  of  the  lemma. 

If  r  or  r'  contains  unbounded  axes,  we  reduce  the  analysis  to  the  bounded  arcs 
by  drawing  a  sufficiently  large  circle  which  contains  all  bounded  axes  and  intersec- 
tion points  of  r  U  r',  cuts  each  unbounded  axe  in  two  points,  and  traverses  only 
unbounded  faces  of  >l(r),  A{T').  The  "chpped"  arrangement  now  consists  of  only 
bounded  arcs,  and  the  unbounded  face  of  each  of  them  is  the  union  of  all  unbounded 
faces  in  the  original  arrangement.  Applying  the  preceding  analysis,  the  claim  follows 
easily  in  this  case  too. 

If  r  U  r'  has  t  unbounded  arcs,  then  there  are  at  most  2(m  +  n  —  t)  endpoints 
and  the  total  nvunber  of  unbounded  faces  is  at  most  2t.  Therefore,  the  above  lemma 
shows  that  it  is  enough  to  compute  only  2(m  +  n-t)  +  2t  =  2(m  -|-  n)  faces  of  each 
of  A{T)  and  A{T')  instead  of  computing  the  whole  arrangement  which  can  have 
Q,{'n?)  (or  fi(m^)  )  faces.  Let  I3{k,n)  (resp.  0'{k,m))  denote  the  maximum  number 
of  edges  in  k  such  faces,  and  let  A{k,n)  (resp.  A'{k,m))  denote  the  time  required 
to  compute  k  faces  of  A{r)  (resp.  A{T')).  We  can  easily  show 

Theorem  4.8  Let  F  be  a  collection  of  n  red  Jordan  arcs  and  let  V  he  a  collection 

of  m  blue  Jordan  arcs.    Then  an  intersection  between  T  and  T'  can  be  detected  in 

time 

O  (  A(2(m  +  n),  n)  +  A'(2(m  -t-  n),  m)+ 

/?(2(m  -I-  n),  n)  -|-  /3'{2{m  +  n),  m))  log(m  -|-  n)) . 

Proof:  By  the  previous  lemma,  it  suffices  to  compute  2{m  +  n)  faces  of  A{T)  and  of 
A{r'),  which  can  be  done  in  time  A(2(m  -|-  n),n)  and  A(2(m  -|-  n),m)  respectively. 
Once  having  computed  these  faces,  an  intersection  between  them  can  be  detected 
in  time  0((/9(2(m  -|-  n),n)  +  ^'{2{m  +  n),m))log(m  -|-  n))  by  using  a  line  sweep. 
Thus,  the  overall  running  time  is  bounded  by 


O  (  A(2(m  -I-  n),  n)  +  A'(2(m  +  n),  m)+ 

0{2{m  +  n),  n)  +  0'{2{m  +  n),  m))  log(m  +  n)) 
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[CEGSW88],  [EGS88],  [EGH*88]  have  studied  the  complexity  of  many  faces  in 
arrangements  of  Jordan  curves  and  for  some  special  ceises  they  have  obtained  sharp 
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(and  sometimes  optimaJ)  bounds,  together  with  comparably  efficient  randomized 
cilgorithms  for  calculating  many  such  faces,  which  are  based  on  the  random  sampling 
technique  of  [C1S7],  [HW87].  Using  these  results  and  Theorem  4.8,  we  obtain 

Corollary  4.9  If  each  ofT,  T'  is  a  set  of  line  segments,  of  unit  circles,  or  of  pseudo 
lines  (i.e.  unbounded  arcs,  each  pair  of  which  intersect  at  most  once),  then  an 
intersection  between  F  and  T'  can  be  detected  in  (randomized)  time  0((m  +  n)''/"'''"'), 
for  any  e  >  0  (where  the  constant  of  proportionality  depends  on  e).  If  T  or  T'  is  a 
set  of  arbitrary  circles,  and  the  other  set  is  either  a  set  of  curves  as  above,  then  an 
intersection  between  T  and  F'  can  be  detected  in  (randomized)  time  0((m  +  n)^'^''"'), 
for  any  e  >  0. 

Proof:  The  proof  follows  immediately  from  the  results  of  [EGS88],  [CEGSW88]  and 
[GSS88].  For  hne  segments,  it  was  shown  in  [EGS88]  that  m  faces  in  an  arrangement 
of  n  line  segments  have  at  most  O  fma'^n^'''^'  +  na{n)j  edges,  and  that  they  can 
be  computed  in  (randomized)  time 

O  ((m3-'n3+2'  +  na(n))  log^  n)  , 
for  any  e  >  0.  In  this  caise  we  have 

A{m  +  n,n)  =  A'(m  +  n,m)  = 

(/3(2(m  +  n),  n)  +  /3'(2(m  +  n),  m))  log(m  +  n)  = 

0{{m  +  n)^/3+«) 

for  any  c  >  0,  so  Theorem  4.8  implies  the  claim.  Similar  bounds  for  arrangements 

of  the  other  kind  of  curves,  as  established  in  [CEGSW88],  [GSS88]  imply  the  claim 

in  cdl  other  cases. 

D 

If  the  axes  in  F  axe  non-intersecting  then  -4(F)  has  only  one  face  and  there  are 
only  n  edges  in  it.  Moreover,  it  is  enough  to  detect  an  intersection  between  F  and 
the  faces  of  A{T')  which  are  unbounded  or  contain  the  endpoints  of  (bounded)  arcs 
in  F.  There  are  at  most  2n  faces  of  ^(F')  containing  endpoints  of  axes  in  F;  the 
unbounded  faces  of  -A(F')  can  be  regcu'ded  as  a  single  connected  face  if  we  clip  each 
unbounded  arc  of  F'  at  two  sufficiently  distant  points,  as  in  the  proof  of  Lemma  4.7. 
We  thus  obtain 

Corollary  4.10  If  the  arcs  in  F  are  non-intersecting,  then  an  intersection  between 
F  and  F'  can  be  detected  in  time  0(A'(2n,m)  +  /3'(2n,m)log(m  +  n)),  where  A'(n,m) 
is  the  time  needed  to  compute  n  faces  of  A{T'),  and  ^'{n,m)  is  the  maximum  number 
of  edges  in  n  such  faces. 
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4.4     Applications 

In  this  subsection,  we  mention  some  applications  of  the  general  red-blue  intersec- 
tion detection  algorithms.    We  describe  only  two  applications  and  leave  it  to  the 
imagination  of  the  reader  to  roam  free  in  seaxch  of  further  applications. 
Consider  the  following  problem: 

Given  a  billiard  table  in  the  shape  of  an  arbitrary  m-sided  convex  poly- 
gon, and  n  billiard  balls  (of  equal  size)  lying  on  the  table.  Suppose  we 
shoot  one  of  these  balls  B  in  a  given  direction,  and  allow  it  to  bounce  k 
times  off  the  sides  of  the  table.   Will  B  hit  any  of  the  other  balls? 

This  problem  can  be  easily  solved  eis  follows.  Expand  each  static  ball  Bi  by  the 
radius  of  B  to  obtain  an  enlarged  ball  B*,  and  shrink  down  the  table  by  the  same 
radius,  so  now  B  becomes  a  moving  point.  The  orbit  of  B  can  easily  be  obtained  in 

time  0{k  logm).  The  results  of  [KLPS86]  imply  that  3  (  Q  B*  j  has  only  0(ti)  arcs 

(note  that  the  expanded  5*'s  can  intersect  each  other).  Now,  we  can  reduce  this 
problem  to  an  instance  of  the  red-blue  intersection  detection  problem  by  considering 

the  orbit  of  B  as  a  collection  T  o{  k  red  circs  and  5  (  [J  B*  j  as  a  collection  T'  of  0{n) 

blue  axes.  Therefore,  it  follows  from  Theorem  4.5  that  we  can  solve  this  problem 
in  time  0{{kJ\^{ri)  +  nJka{k))\og^'^{n  +  k)).  However,  in  this  problem  the  arcs 
satisfy  the  following  additional  properties.  The  orbit  of  B  forms  a  connected  chain 
of  line  segments,  i.e.  F  can  be  partitioned  into  subsets  of  airbitrary  size  so  that  the 
boundary  of  each  face  in  the  arrangement  of  each  subset  is  connected.  Also,  the  arcs 
in  r'  are  non-intersecting.  Hence,  it  follows  from  the  remarks  following  Theorem  4.5 
that  we  can  solve  this  problem  in  time  0{kJX4{n)  log*'^  n)  or  in  0{n\/k\ogn).  Our 
solution  can  easily  be  generalized  to  hzmdle  billiard  "balls"  of  more  general  shape. 
These  properties  also  imply  that  we  can  use  Corollary  4.10  to  check  whether  B  hits 
any  static  ball  in  (randomized)  time  0((A:-|-n)''/^''"'),  for  any  e  >  0.  (The  assumption 
that  the  table  is  convex  was  made  to  ensure  that  the  shrunk  table  is  still  polygonal, 
it  can  probably  be  relaxed  to  allow  any  simple  polygonal  table,  although  in  this 
case  the  shrunk  table  may  have  some  circular  eltcs  on  its  boundary,  in  which  case 
some  adaption  of  the  ray  shooting  technique  of  [CG85]  will  be  required  to  obtain 
efficiently  the  orbit  of  B.) 

Another  application,  more  related  to  computer  graphics,  is  multiple  ray-shooting 
(the  "death  squad"  problem): 

Given  m  rays  in  the  plane  and  a  collection  of  objects  bounded  by  n 
algebraSc  arcs  (of  small  degree),  does  any  ray  hit  any  object  ? 
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Again  this  problem  caji  be  reduced  to  an  instance  of  the  red-blue  intersection  de- 
tection problem  by  considering  the  boundaries  of  the  objects  as  a  collection  of  red 
arcs  and  the  rays  as  a  collection  of  blue  arcs.  Hence,  assuming  the  objects  axe 
non-intersecting,  we  can  solve  this  problem  in  (deterministic)  time  0(m-y/n  log  n) 
(by  Remark  (ii)  following  Theorem  4.5),  or  in  randomized  time  0(n^/^~^'m^/^''"'  -|- 
mQ(m)log^  m),  for  any  e  >  0  (using  Corollary  4.10). 

If  instead  of  rays  we  shoot  (identiced)  bullets  of  some  shape,  we  need  to  expand 
each  object  by  the  bullet  shape  to  reduce  the  problem  to  the  original  ray  shooting 
problem.  In  this  case  we  can  detect  em  intersection,  using  our  general  algorithm,  in 
time  0({rnJXg+2in)  +  nJrna{m))log^'^{m  +  n)),  where  s  is  the  maximum  number  of 
intersections  between  amy  two  expanded  objects.  This  can  be  improved  by  noting 
that  a  single  face  in  an  arrangement  of  m  rays  has  only  0{m)  complexity  (and 
can  be  calculated  in  O(mlog^m)  time).  Moreover,  if  objects  and  the  bullets  are 
all  convex,  then  s  =  2  [KLPS86],  and,  talcing  expanded  object  as  a  closed  Jordan 
curve,  we  obtEiin  a  0{(my/n  +  ny/Tn)\og^'^{m  +  n))  solution. 

Our  techniques  also  allow  us  to  solve  efficiently  various  other  extensions  of  the 
above  problem.  For  example,  if  the  bullets  being  shot  have  a  finite  range,  then 
we  need  to  detect  axi  intersection  between  a  collection  of  segments  and  the  given 
objects.  Similarly,  in  the  case  of  a  "stone-throwing  squad",  where  all  the  stones 
move  along  parabolic  trajectories  in  the,  say  x  —  z  plane,  or  along  other  (algebraic) 
trajectories,  we  need  to  detect  an  intersection  between  those  trajectories  emd  the 
objects. 

5      Conclusions 

In  this  paper  we  obtained  severed  efficient  algorithms  for  the  red-blue  intersection 
detection  problem.  Our  algorithm  for  the  case  when  each  of  A{T)  and  A{T')  is 
a  connected  planar  graph  is  close  to  optimal,  and  our  general  (deterministic  and 
randomized)  algorithms  are  faster  than  previously  known  algorithms.  We  applied 
these  algorithms  to  obtain  fast  algorithms  for  several  other  problems  in  motion 
planning  and  collision  detection.  However,  there  axe  still  several  open  questions. 
The  most  importemt  is  whether  the  running  time  of  our  deterministic  algorithm  for 
the  general  case  can  be  improved.  Theorem  4.8  shows  that  an  efficient  algorithm  for 
computing  many  faces  in  the  arrangement  of  a  given  set  of  Jordan  arcs  also  yields  a 
fast  algorithm  for  the  red-blue  intersection  detection  problem,  therefore  one  way  of 
improving  the  running  time  is  to  find  an  efficient  deterministic  algorithm  to  compute 
many  faces  in  the  arremgement  of  Jordan  axes.  But  this  problem  in  general  seems 
to  be  very  haxd,  so  it  would  be  interesting  to  obtain  efficient  cdgorithms  even  for 
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special  cases  like  line  segments,  etc. 

We  also  obtained  randomized  algorithms  for  several  special  cases,  which  are 
faster  than  our  deterministic  algorithm  but  at  present  we  do  not  have  any  such 
algorithm  for  the  general  case.  Another  open  question  is  whether  there  exist 
0(n  log*^'^^  n)  randomized  algorithms  if  all  the  eircs  are  line  segments.  Finally,  can 
one  report  all  red-blue  intersections  in  an  output-sensitive  fashion  for  general  arcs 
also? 
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Figure  1:  Red  arcs  and  blue  polygon 


Figure  2:  A{T)  and  A{T')  are  planar  networks. 
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Figure  3:  Simple  polygon  inside  a  corridor 


Figure  4a:  A  corner  with  6,  a:  6, 


Figure  4b:  A  corner  with  6,  <  6, 
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Figure  5:  Lemma  4.1 


Figure  6:  Lemma  4.7 
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